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Abstract 

The paradox of destabilization of a conservative or non-conservative system by small dissipation, or 
Ziegler's paradox (1952), has stimulated an ever growing interest in the sensitivity of reversible and Hamil- 
tonian systems with respect to dissipative perturbations. Since the last decade it has been widely accepted 
that dissipation-induced instabilities are closely related to singularities arising on the stability boundary. 
What is less known is that the first complete explanation of Ziegler's paradox by means of the Whitney 
umbrella singularity dates back to 1956. We revisit this undeservedly forgotten pioneering result by Oene 
Bottema that outstripped later findings for about half a century. We discuss subsequent developments of 
the perturbation analysis of dissipation-induced instabilities and applications over this period, involving 
structural stability of matrices, Krein collision, Hamilton-Hopf bifurcation and related bifurcations. 
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1 Introduction 

'II n'y a de nouveau que ce qui est oublie' — this paraphrase of the Ecclesiastes 1:10, attributed to Marie- 
Antoinette, perfectly summarizes the story of the mathematical description of the destabihzing effect of vanishing 
dissipation in non-conservative systems. 
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There is a fascinating category of mechanical and physical systems which exhibit the following paradoxical 
behavior: when modeled as systems without damping they possess stable equilibria or stable steady motions, 
but when small damping is introduced, some of these equilibria or steady motions become unstable. 

The paradoxical effect of damping on dynamic instability was noticed first for rotor systems which have 
stable steady motions for a certain range of speed, but which become unstable when the speed is changed to a 
value outside the range. 

In 1879 Thomson and Tait \7W showed that a statically unstable conservative system which has been sta- 
bilized by gyroscopic forces could be destabilized again by the introduction of small damping forces. More 
generally, they consider conservative and nonconservative linear two degrees of freedom systems in remarkable 
detail. The destabilization by damping, using Routh's theorems, is implicit in their calculations, it is not 
formulated as a paradox. 

In 1924, to explain the destabilization of a flexible rotor in stable rotation at a speed above the critical 
speed for resonance, Kimball |36| introduced a damping of the rotation, which has lead to non-conservative 
positional (circulatory) forces in the equations of motion of a gyroscopic system. In 1933 Smith [76] found that 
this non-conservative rotor system loses stability when the speed of rotation ft > w(l + -^), where uj is the 
undamped natural whirling frequency (the critical speed for resonance) and 6 and v are the viscous damping 
constants for the stationary and rotating damping mechanisms. In Smith's model, the destabilizing effect of 
the damping of rotation (v), observed also by Kapitsa [34], was compensated by the stationary damping (5). 
This was a first demonstration of a strong influence of the spatial distribution of damping (or equivalently the 
modal distribution) on the borderlines between stability and instability domains in multi- modal non-conservative 
systems [l9 l [89 ] . 

In the 1950s and 1960s the publications of Ziegler [90l[91], Bolotin [IHlin], and Herrmann [25l[26], motivated 
by aerodynamics applications, initiated a considerable activity in the investigation of dynamic instability of 
equilibrium configurations of structures under non-conservative loads. The canonical problem was the flutter 
of a vertical flexible cantilever column under a compressive non-conservative or follower load which remains 
tangent to the bending column. In the flutter mode the tip of the column is preponderantly slanted towards the 
left during the half-cycle in which the tip is moving towards the right and vice versa in the following half-cycle. 
This snake-like oscillation permits the follower force to do positive work on each cycle ^T^ . 

The strong influence of the spatial or modal distribution of damping within the structure on its stability 
under non-conservative loading, observed in these publications, should not have been surprising in the light 
of earlier findings of rotor dynamists. However, they revealed explicitly the most dramatic and paradoxical 
aspect of the sensitivity of the stability of the nonconservative structures to small damping forces. It turned out 
that the critical load for a structure with small damping may be considerably smaller than that for the same 
structure without damping. In other words, there is a wide range of loads for which the undamped structure is 
stable, but which produce instability as soon as a tiny bit of damping is added to the structure. 

These phenomena were actively studied in the 1960's to provide more basic understanding and they have 
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continued to be studied with more sophisticated tools, including early attempts to employ singularity theory 
[5T] . until in the mid 1990s it was understood [SUIIZI] that the destabilization paradox is related to the Whitney 
umbrella singularity of the stability boundary [MUSS]. After describing in the first sections Whitney's umbrella 
and Ziegler's paradox, we make in scction[4]a sharp turn to the 1950s to revisit an article of Oene Bottema jT4] . 
who in 1956 first made this discovery and clarified the paradox. Surprisingly, this paper surpassed the attention 
of most scientists during five decades. 

In section [2] we will relate these results to singularity theory, in sections [5] and [7] we show in various ways 
their extension to finite- and infinite-dimensional systems using perturbation theory of multiple eigenvalues, in 
section[B]we focus on periodic systems, and in the remainder we discuss applications in physics and engineering. 

2 Whitney's umbrella 

In a remarkable paper of 1943 [84], Hassler Whitney described singularities of maps of a differential n-manifold 
into with m = 2n — 1. It turns out that in this case a special kind of singularity plays a prominent 
role. Later, the local geometric structure of the manifold near the singularity has been aptly called 'Whitney's 
umbrella'. In Fig. [T] we reproduce the original sketch of the singular surface from the companion article [55] . 
The paper contains two main theorems. Consider the C'^ map / : E"^ with m = 2n — 1. 

1. The map / can be altered slightly, forming /*, for which the singular points are isolated. For each such 
an isolated singular point p, a technical regularity condition C is valid which relates to the map /* of the 
independent vectors near p and of the differentials, the vectors in tangent space. 

2. Consider the map /* which satisfies condition C. Then we can choose coordinates x — {xi, X2, - ■ ■ , Xn) in 
a neighborhood of p and coordinates y = {yi, 2/2, ■• ■ , ym) (with m 2n — 1) in a neighborhood of y = /(p) 
such that in a neighborhood of f*{p) we have exactly 

2/1 = xl^ 

yi = Xi, i = 2, - • • ,n, 
Vn+i-i = xiXi, j = 2, •••,n. 

If for instance n = 2, m = 3, the simplest interesting case, we have near the origin 

2/1 = 2/2 = X2, 2/3 = XiX2, (1) 

so that 2/1^0 and on eliminating xi and X2'- 

2/12/1 - 2/3 = 0- (2) 

Starting on the 2/2-a'Xis for 2/1 = 2/3 = 0, the surface widens up for increasing values of 2/1- For each 2/2, the 
cross-section is a parabola; as 2/2 passes through 0, the parabola degenerates to a half-ray, and opens out again 
(with sense reversed); see Fig.[T] 
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Figure 1: Whitney's original sketch of the umbrella [85] , 

Note that because of the C'^ assumption for the differentiable map /, the analysis is delicate. There is a 
considerable simplification of the treatment if the map is analytical. 

The analysis of singularities of functions and maps is a fundamental ingredient for bifurcation studies of 
differential equations. After the pioneering work of Hassler Whitney and Marston Morse, it has become a huge 
research field, both in theoretical investigations and in applications. We can not even present a summary of 
this field here, so we restrict ourselves to citing a number of survey texts and discussing a few key concepts and 
examples. In particular we mention [5], [5T], [55], [J], [5] and [S]. A monograph relating bifurcation theory with 
normal forms and numerics is [55] . 

The relation between singularities of functions and critical points or equilibria of differential equations 
becomes relatively simple when considering Hamiltonian and gradient systems. Consider for instance the time- 
independent Hamiltonian function H{p, q) with p, q G K". Singularities of the function H are found in the set 
R2" where 

dH _dH _^ 
dp dq 

These points correspond with the critical points (equilibria) of the Hamiltonian equations of motion 

. dH . _ dH 
^ 9q ' ^ dp 

More in general, consider the dynamical system described by the autonomous ODE 

x = f(x), xgR", f :R"h^R". 
An equilibrium xq of the system arises if f (xq) — 0. With a little smoothness of the map f we can linearize 
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near Xq so that we can write 



X 



A(x - xo) + g(x) 



(3) 



with A a constant n x n— matrix, g(x) contains higher-order terms only. In other words 




g(x) is tangent to the Unear map in xq. 

The properties of the matrix A determine in a large number of cases the local behavior of the dynamical 
system. In a seminal paper [3 , Arnold considers families of matrices, smoothly depending on a number of 
parameters (denoted by vector p). So, for the constant n x n— matrix we write Ap. Suppose that for p = 0, Aq 
is in Jordan normal form. Choosing p in a neighborhood of p = produces a deformation (or perturbation) of 
Aq, assuming that near p = the entries of Ap can be expanded in a convergent power series in the parameters. 
A deformation is versal if all other deformations near p = are equivalent under smooth change of parameters. 

The paper [3^ uses normal forms to obtain suitable versal deformations. These are associated with the 
bifurcations of the linearized system ([3]). Note that although a matrix induces a linear map, the corresponding 
eigenvalue problem produces a nonlinear characteristic equation. In addition, the parameters involved, make 
it necessary to analyze maps of R" into M™. For instance in the following sections we meet with maps from 
M? into R-^ as studied by Whitney [84]. Nevertheless, in 1943 it was hard to imagine that this study of global 
analysis, a pure mathematical abstraction, would find already an industrial application in the next decade. 

3 Ziegler's paradox 

In 1952 Hans Ziegler of ETH Zurich published a paper [90] that became classical and widely known in the 
community of mechanical engineers; it also attracted the attention of mathematicians. Studying a simplified two- 
dimensional model of an elastic rod, fixed at one end and compressed by a tangential end load, he unexpectedly 
encountered a phenomenon with a paradoxal character: the domain of stability of the Ziegler's pendulum 
changes in a discontinuous way when one passes from the case in which the damping is very small to that where 
it has vanished [501 [OT] . 

Ziegler's double pendulum presented in Fig.[2Ka) consists of two rigid rods of length I each, whose inclinations 
with respect to the vertical are denoted as Lpi and ^p2- Two masses mi and m2 with the weights Gi and G2 are 
concentrated at the distances oi and 02 from the joints. The elastic restoring torques and the damping torques 
at the joints are api, c{(p2 — ^1) and biLpi, b2{'f2 — 'Pi)j respectively. With these assumptions the kinetic energy 
of the system is 




[{■miol + ■m2l^)ipi + 2m2la2ipiip2 + 111-20^(^1] , 



(4) 



while the potential energy reads 



^ = 2 [(Giai + G2I + 2c)lpI - 2c(^i(^2 + (6202 + 0)^32] ■ 



(5) 
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a) 




o(Vs-9i}*b(h-'Pi) 




Figure 2: Original drawings from the Ziegier's work of 1952 iDO]: (a) double linked pendulum under the follower 
load, (b) (bold line) stability interval of the undamped pendulum and (shaded area) the domain of asymptotic 
stability of the damped pendulum with equal coefficients of dissipation in both joints. If 6 = we have no 
dissipation and stability for follower force P \i P < Pk. 



The generalized dissipative and non-conservative forces are then 

Qi = Pl(ip\ - ip2) - iibi + b2)ipi ~ 62(^2), Q2 = b2{(pi - (P2)- 



(6) 



Writing the Lagrange's equations of motion 



Qi, where L = T — V and a dot denotes time 



differentiation, and assuming Gi = and G2 — for simplicity, we find 



miol + m2p TO2?a2 



7712^02 17120,2 





= 0. 



(7) 



With the substitution (pi ~ Aiexp{\t), equation ([7]) yields a 4-dimensional eigenvalue problem with respect to 
the spectral parameter A. 

Putting mi = 2m, m2 = to, ai 02 = ^, &i = 62 = 6 and assuming that dissipation is absent (6 ~ 0), 
Ziegler found from the characteristic equation that the vertical equilibrium position of the pendulum looses its 
stability when the magnitude of the follower force exceeds the critical value Pk , where 

'7 



Pk 



--^/2^. 2.086f 



(8) 



In the presence of damping {b > 0) the Routh-Hurwitz condition yields the new critical follower load that 
depends on the square of the damping coefficient b 



/,x 41c 1 52 
^^^(')= 287 + 2;^- 



(9) 



Ziegler found that the domain of asymptotic stability for the damped pendulum is given by the inequalities 
P < Pk{b) and b > and he plotted it against the stability interval P < Pk of the undamped system. Fig. [21(b). 
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Surprisingly, the limit of the critical load Pk {b) when b tends to zero turned out to be significantly lower than 
the critical load of the undamped system 



Note that in the original work of Ziegler, formula ([9]) contains a misprint which yields linear dependency of 
the critical follower load on the damping coefficient b. Nevertheless, the domain of asymptotic stability found 
in [SD] and reproduced in Fig. O^b), is correct. 

Some authors considered extensions of Ziegler's model by adding a conservative load and by assuming unequal 
damping coefficients [11] [351 [33 [SD [HH] ■ Fig. [3] demonstrates how the domain of instability for the undamped 
Ziegler's pendulum with the partially follower load (rj — 1 corresponds to the pure follower load), shown in dark 
gray in the (?7,p)-plane, extends in a discontinuous manner in the presence of dissipation when 62 = O.Sfei and 
bi 0. The portion of the stability domain that became unstable is depicted in light gray [37l[80]. Therefore, 
the two-dimensional stability diagrams of the undamped system and the system with vanishingly small damping 
differ by a region of positive measure. 

Ziegler drew attention both to the substantial decrease in the critical load of the damped non-conservative 
system with vanishingly small dissipation and to the high sensitivity of the critical follower load with respect to 
the variation of the damping distribution. In the mechanical engineering literature these two effects are called 
the Ziegler's paradox of destabilization by small damping. 



Figure 3: Ziegler's pendulum with the partially follower force [37]: (dark gray) Instability domain in the absence 
of damping (61 = 0, 62 = 0) and (light gray) its increment in the presence of damping with &2 — 0.36i and bi — s- 0. 

In the conclusion to his classical book [TU], Bolotin emphasized that the discrepancy between the stability 
domains of undamped non-conservative systems and that of systems with infinitesimally small dissipation is a 
topic of the greatest theoretical interest in stability theory. Encouraging further research of the destabilization 
paradox, Bolotin was not aware that the crucial ideas for its explanation were formulated as early as 1956. 



(10) 
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4 Bottema's solution 



In 1956, in the journal 'Indagationes Mathematicae', there appeared an article by Oene Bottema (1901-1992) 
[H] , then Rector Magnificus of the Technical University of Delft and an expert in classical geometry and me- 
chanics, that outstripped later findings for decades. Bottema's work in 1955 [13 can be seen as an introduction, 
it was directly motivated by Ziegler's paradox. However, instead of examining the particular model of Ziegler, 
he studied in [H] a much more general class of non-conservative systems. 

Following |131 114j , we consider a holonomic scleronomic linear system with two degrees of freedom, of which 
the coordinates x and y are chosen in such a way that the kinetic energy is T = (i^ + j/^)/2. Hence the Lagrange 
equations of small oscillations near the equilibrium configuration x ~ y — are as follows 

X + aiix + ai2y + biix + bi2y = 0, 

y + a2ix + a22y + b2ix + b22y = 0, (11) 

where and bij are constants, A := (oy ) is the matrix of the forces depending on the coordinates, B :— {btj) 
of those depending on the velocities. If A is symmetrical and while disregarding the damping associated with 
the matrix B, there exists a potential energy function V — (anx^ -I- 2ai2xy + a22y^) /2, if it is antisymmetrical, 
the forces are circulatory. When the matrix B is symmetrical, we have a non-gyroscopic damping force, which 
is positive when the dissipative function (feux^ -I- 2bi2xy + 622y^)/2 is positive definite. If B is antisymmetrical 
the forces depending on the velocities are purely gyroscopic. 

The matrices A and B can both be written uniquely as the sum of symmetrical and antisymmetrical parts: 
A = K + N and B = D + G, where 

/ fell fci2 \ f v\ / dii di2 \ f n\ 

K = " , N = , D = " , G = , (12) 

\ fc21 ^22/ J \ d21 d22 ) \-n J 

with fell = On, fc22 = 022, ki2 = k2i = (ai2 + a2i)/2, v = (ai2 - a2i)/2 and dn = bn, ^22 = ^^22, ^12 = ^21 = 

(&12+&2l)/2, f^- (fel2-fo2l)/2. 

The system (|12|) has a potential energy function (disregarding damping) when = 0, it is purely circulatory 
for fell — ki2 = ^22 = 0, it is non-gyroscopic for il = 0, and has no damping when dn = di2 = (^22 = 0. If 
damping exists, we suppose in this section that it is positive. 

In order to solve the equations (I12p we put x ~ Ciexp(Ai), y ~ C2exp(Ai) and obtain the characteristic 
equation for the frequencies of the small oscillations around equilibrium 

Q A'' + aiA^ + a2A2 +a3A + a4 = 0, (13) 

where [SHISHISS] 

ai=trD, a2 = trK + det D + 03 = trKtrD - trKD + 2f7i^, a4 = dctK + iy2. (14) 

For the equilibrium to be stable all roots of the characteristic equation (jl3p must be semi-simple and have real 
parts which are non positive. 
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It is always possible to write, in at least one way, the left hand-side as the product of two quadratic forms 
with real coefficients, Q = (A^ +PiA + (7i)(A^ +^2^^ + 92)- Hence 

ai=Pi+P2, 0,2 ^ P1P2 + qi + q2, 03 =Pi92 +^291, ai=qiq2- (15) 

For all the roots of the equation (fT3|) to be in the left side of the complex plane {L) it is obviously necessary 
and sufficient that pi and qi are positive. Therefore in view of (llSp we have: a necessary condition for the roots 
Q = having negative real parts is ai > {i — 1,2,3,4). This system of conditions however is not sufficient, 
as the example (A^ - A + 2){\^ + 2A + 3) = A'' + A^ + 3A^ + A + 6 shows. But if > it is not possible 
that either one root of three roots lies in L (for then a4 < 0); it is also impossible that no root is in it (for, 
then 04 < 0). Hence if > at least two roots are in i; the other ones are either both in L, or both on the 
imaginary axis, or both in R. In order to distinguish between these cases we deduce the condition for two roots 
being on the imaginary axis. If p.i {p. ^ is real) is a root, then /i"' — a2/i^ + 04 = and — ai/i^ + 03 = 0. Hence 
H := altti + flg — 010203 = 0. Now by means of psp we have 

H = -piP2{aia3 + (qi ~ q2)^)- (16) 

In view of oi > 0, 03 > the second factor is positive; furthermore oi = pi +P2 > 0, hence pi and p2 cannot 
both be negative. Therefore H < implies pi > 0, p2 > 0, for H ~ we have either pi — or P2 = (and 
not both, because 03 > 0), for iJ > pi and p2 have different signs. We see from the decomposition of the 
polynomial (jl3p that all its roots are in L if pi and p2 are positive. 

Hence: a set of necessary and sufficient conditions for all roots of (fTS]) to be on the left hand-side of the 
complex plane is 

a,; >0 {i = 1,2,3,4), H <0. (17) 

We now proceed to the cases where all roots have non-positive real parts, so that they lie either in L or on 
the imaginary axis. If three roots are in L and one on the imaginary axis, this root must be A = 0. Reasoning 
along the same lines as before we find that necessary and sufficient conditions for this are > {i — 1,2,3), 
04 = 0, and H < 0. If two roots are in L and two (different) roots on the imaginary axis we have pi > 0, qi > 0, 
P2 = 0, q2 > and the conditions are > (i = 1, 2, 3, 4) and H = 0. If one root is in L and three are on 
the imaginary axis, then pi > 0, qi — 0, P2 = 0, 92 > and the conditions are > {i — 1, 2, 3), 04 — 0, and 
H = 0. 

The obtained conditions are border cases of (flT)) . This does not occur with the last type we have to consider: 
all roots are on the imaginary axis. We now have pi = 0, P2 = 0, qi > 0, q2 > 0. Hence 02 > 0, 04 > 0, 
ai — as ~ and therefore H — 0. This set of relations is necessary, but not sufficient, as the example 
Q — + 6A^ + 25 = (which has two roots in L and two in the righthand side of the complex plane (i?)) 
shows. The proof given above is not valid because as seen from (jl7p . H — does not imply now piP2 = 0, the 
second factor being zero for aia^ — and qi — q2- The condition can of course easily be given; the equation 
(fT3| is A^ + a2)? -f 04 = and therefore it reads 02 > 0, 04 > 0, a\ > Aa^. 
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Summing up we have: all roots of (|13p (assumed to be different) have non-positive real parts if and only if 
one of the two following sets of conditions is satisfied [13] 

a?a4 + a? 

A : ai > 0, 02 > 0, 03 > 0, 04 > 0, a2 > — -, 

B : fli = 0, 02 > 0, 03 = 0, 04 > 0, a2 > 2^/04. (18) 

Note that ai represents the damping coefficients 611 and 622 in the system. One could expect _B to be a 
limit of A, so that for oi ^ 0, 03 ^ the set A would continuously tend to B. That is not the case. 

Remark first of all that the roots of (fT3|) never lie outside i? if ai = 0, 03 7^ (or ai 7^ 0, 03 ~ 0). 
Furthermore, if A is satisfied and we take ai — e6i, 03 — sbs, where bi and 63 are fixed and £ ^ 0, the last 
condition of A reads (e ^ 0) 

6?a4 + 63 

«2 > , , = 91 
bibs 



while for e = we have 



Obviously we have [Mj 



02 > 2\/aI = 32- 
(foiVoI- b 

31 - 52 = ■ 



bibs 

so that {gi > 172) but for 63 — biy/a^. That means that in all cases where 63 ^ biy/a^ we have a discontinuity 
in our stability condition. The phenomenon of the discontinuity was illustrated by Bottema in a geometrical 
diagram, Fig. H) 

Following Bottema we substitute in (fTH]) A ~ cfi, where c is the positive fourth root of 04 > 0. The new 
equation reads P := /i^ + fei/i'' + 62^^ + bsfJ. + 1 = 0, where bi = ai/d {i — 1, 2, 3, 4). If we substitute — c^bi in 
A and B we get the same condition as when we write bi for a^, which was to be expected, because if the roots 
of (|13p are outside R, those of P = are also outside R and inversely. We can therefore restrict ourselves to the 
case 04 = 1, so that we have only three parameters oi, 02, 03. We take them as coordinates in an orthogonal 
coordinate system. 

The condition H = 01 

010203 = 0^+03 (19) 

is the equation of a surface V of the third degree, which we have to consider for oi > 0, 03 > 0, Fig. U) Obviously 
1^ is a ruled surface, the line 03 = moi, 02 = m + 1/ni (0 < m < 00) being on V. The line is parallel to the 
0oiO3-plane and intersects the 02-axis in oi = 03 = 0, 02 = m + 1/m > 2. The 02-axis is the double line of V, 
02 > 2 being its active part. Two generators pass through each point of it; they coincide for 02 = 2 (m = 1), 
and for 02 — » 00 their directions tend to those of the oi and a3-axis (m = 0, to = 00). The conditions A and B 
express that the image point (oi, 02, 03) lies on V or above V. The point (0, 2, 0) is on V, but if we go to the 
a2-axis along the line 03 = tooi the coordinate 02 has the limit m+ 1/m, which is > 2 but for m = 1. Curiously 
enough, even half a century later, there appear papers repeating this reasoning and the result almost literally, 
see for instance 1551. 
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Figure 4: Original drawing (left) from the 1956 work 14J of Oene Bottema (right), showing the domain of the 
asymptotic stability of the real polynomial of fourth order and of the two-dimensional non-conservative system 
with Whitney's umbrella singularity. The ruled surface (called V in the text) is given by equation ([TI 



Note that we started off with 8 parameters in eq. (fTT|) , but that the surface V bounding the stability domain 
is described by 3 parameters. It is described by a map of E'^ into E'^ as in Whitney's papers [SHIHS]. Explicitly, 
a transformation of (19) to (2) is given by 

1 „ 1 

ai = -ys + w, a2 = 2 + y2, ~-y3 + w 

Wfth = +y^y^_ 

Returning to the non-conservative system pT|) {ly ^ 0), with damping, but without gyroscopic forces, so 
57 = 0, and assuming as in 13J that ki2 = 0, fcn > 0, and fc22 > (a similar setting but with di2 and 
ki2 7^ was considered later by Bolotin in [10l[T^), we find that the condition for stability H <0 reads 

2 (fell - fc22)^ _ {dll - ^22)^(^11 - fc22)^ - 4(^11^22 + ^22 C^ll ) (c^ll ^^22 " C?12)(^ll + ^22) .„„,^ 

^ 4 4(dn+d22)2 ■ ^ ' 

Suppose now that the damping force decreases in a uniform way, so we put dn = ed[-^, di2 — sd'12, ^22 — £^221 
where dn, (ii2, ^22 are constants and e — > 0. Then, for the inequality (pO)) we get 



2^ 2 _ (fell ~ fc22)^ «! - d'22?{kii - k22? 



4 mi+d'22? 
But if there is no damping, we have to make use of condition B, which gives 



o n (fell - fc22)^ ArK\~ , , , 

<^p2 :^ ^ " ^ ^ f — j -detK. (22) 



11 



Obviously 



2 „2 _ id'll - ^22)^(^11 - "-22; _ z,u,xvx^ - uxxvuxx^ 

where the expressions written in terms of the invariants of the matrices involved [45] are valid also without the 
restrictions on the matrices D and K that were adopted in [lOl [TH [13] . For the values of which 
are small with respect to uq we can approximately write 30] 



2)' 



2trKD trKtrD 



1^0 



1 
2^0 



2trKD - trKtrD 



2trD 



(24) 



If D depends on two parameters, say Si and S2, then ([M)) has a canonical form ([2]) for the Whitney's umbrella 
in the ((5i, (52, i^)-space. Due to discontinuity existing for 2trKD — trKtrD 7^ the equilibrium may be stable 
if there is no damping, but unstable if there is damping, however small it may be. We observe also that the 
critical non-conservative parameter, i/cr, depends on the ratio of the damping coefficients and thus is strongly 
sensitive to the distribution of damping similarly to how it happens in rotor dynamics. This is the results which 
Ziegler [5D1 [OT] found in a special case. 



5 'Hopf meets Hamilton under Whitney's umbrella' 

The title of this section derives from a nice tutorial paper by Langford |56j . As we have seen, Bottema was 
the first who established that the asymptotic stability domain of a real polynomial of fourth order in the space 
of its coefficients consists of one of the 'pockets' of the Whitney umbrella. The corresponding singularity was 
later identified as generic in the three parameter families of real matrices by V.I. Arnold [31 [4], who named it 
'deadlock of an edge'. In this respect Bottema's results in |14] can be seen as an early study of bifurcations and 
structural stability of polynomials and matrices, and therefore of the singularities of their stability boundaries 
whose systematical treatment was initiated since the beginning of the 1970s in [3l [4l [57l I58j . 

Although Bottema applied his result to nonconservative systems without gyroscopic forces, there are reasons 
for the singularity to appear in the case when gyroscopic forces are taken into account because the stability 
is determined by the roots of a similar fourth order characteristic polynomial. In order to study this case we 
consider separately the following m-dimcnsional version of the non-conservative system (llip 

x+{nG + ST>)±+{K + iyN)x = 0, (25) 

where a dot stands for time differentiation, x e M™, and real matrix K = K-'" corresponds to potential forces. 
Real matrices D = D^, G = -G^, and N = -N^ are related to dissipative (damping), gyroscopic, and non- 
conservative positional (circulatory) forces with magnitudes controlled by scaling factors 6, f2, and ly respectively. 
A circulatory system, to which the undamped Ziegler's pendulum is attributed [401 1621 [70] . is obtained from 
([^5)1 by neglecting velocity-dependent forces 

X + (K + iyN)x = 0, (26) 
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while a gyroscopic one has no damping and non-conservative positional forces 



X + nG± + Kx = 0. 



(27) 



Circulatory and gyroscopic systems (pS)) and (^7)) possess fundamental symmetries that are evident after trans- 
formation of equation (|25p to the form y = Cy with 



C = 



-il7G I 

i^rJDG + ir^^G^ - K - z^N -^D-irJG 



y = 



x+if^Gx 



(28) 



where I is the identity matrix. 

In the absence of damping and gyroscopic forces (5 = = 0), RCR = — C with 



R = R 



I 

-I 



(29) 



This means that the matrix C has a time reversal symmetry, and equation ()26p describes a reversible dynamical 
system |62j . Due to this property, 



det(C - AI) = det(R(C - AI)R) = det(C + AI), 



(30) 



and the eigenvalues of circulatory system (l26l) appear in pairs (—A, A). Without damping and non-conservative 
positional forces {5 = v = Q) the matrix C possesses the Hamiltonian symmetry JCJ = C"^, where J is a 
symplectic matrix [H [501 [H] with 

I 



As a consequence, 



-I 



det(C - AI) = det(J(C - AI)J) = det(C^ + AI) = det(C + AI), 



(31) 



(32) 



which implies that if A is an eigenvalue of C then so is —A, similar to the reversible case. Therefore, an equilibrium 
of a circulatory or of a gyroscopic system is either unstable or all its eigenvalues lie on the imaginary axis of 
the complex plane, in the last case implying marginal stability if they are semi-simple. 

It is well known that in the Hamiltonian case, the transition from gyroscopic stability to flutter instability 
occurs through the interaction of simple purely imaginary eigenvalues with the opposite Krein signature known 
as the Krein collision or the Hamiltonian Hopf bifurcation [23l [56l [59l [60l [61] . The collision occurs at the border 
of marginal stability, say at = fJo for (|27p . and it yields a double pure imaginary eigenvalue with the Jordan 
chain of vectors, which splits into a a complex conjugate pair under destabilizing variation of the parameter fJ. 

Let iujQ be the double eigenvalue at = fJo with the Jordan chain of generalized eigenvectors Uq, Ui, 
satisfying the equations [46] 



{-IijjI + iuqVIqG + K)uo 



0, 



(— Iwq + io^o^^oG + K)ui = — (2ia;oI + ^loG)uo. 



(33) 
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Then, the Krein coUision in the gyroscopic system (|27p is described by the foUowing expressions 

u{n) = UQ±i^lUl^/n~Ck) + o{\^l-no\^), (34) 

where the real coefficient fu, is according to [15] 



_i_ Q rxi rri _i_ . V / 

ilp (ti'oUi ui — u;[ Kui — icjoi^o u;[ Gui — Uq uq) 



with the bar over a symbol denoting complex conjugate. 

Perturbing the system (|27p by small damping and circulatory forces yields an increment to a simple pure 
imaginary eigenvalue [101 US] 

^ = -Ttr,\x^ /n^ , 2-Tm^ ^oa +o{S,i^). (36) 

u (i2)Ku(i2) + cj^u (i2)u(S2) 

With the expressions , equation (|36p is used for the calculation of the deviation from the imaginary axis 
of the eigenvalues that participated in the Krein collision in the presence of the non-Hamiltonian perturbation 
that makes the merging of modes an imperfect one |28j . 

Since D and K are real symmetric matrices and N is a real skew-symmetric one, the first-order increment 
to the eigenvalue iwifl) given by (j36p is real-valued. Consequently, in the first approximation in S and i^, the 
simple eigenvalue iuj{^l) remains on the imaginary axis, ii v = j(fl)6, where 

'u^(J7)Nu(J7)' 

With the expansions ([34|l the formula ((37|l reads 



= ~^u^m :l^'^';^z Z - (37) 



(38) 



where we define 



di — Re(ug Duq), d2 = Im(ug Dui — Duq), 

ni = Im(uoNuo), ^2 = Re(uo Nui - Nuq), = -i^o "r!^"" ■ (39) 

Uq Nuo 

From (|38p it follows that in the vicinity of 7 := = 7, the limit of the critical value of the gyroscopic parameter 
flcr of the near-Hamiltonian system as (5 — s- exceeds the threshold of gyroscopic stabilization determined by 
the Krein collision (see [15]) 

n,rii) = no+ > no. (40) 

Substituting 7 = //(5 in expression (|40p yields a simple estimate for the critical value of the gyroscopic 
parameter ilcr{S, v) that has a canonical form ^ and therefore describes the Whitney's umbrella surface in the 
((5, V, r2)-space 06] 

i7,,,(<5, .) = + ,2» - (41) 

fi-'(uJod2 - 7*n2 - di)^()'^ 
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In case of two oscillators (m = 2) the approximation (j4ip is transformed to 

2 ..2 trKD 



^ 2 , ^,3 trKD + (f^o^ - t^o)*!'^ 

ri,,((5,j.) -r!o + f^o^— F:t7T7(i'-7*'5)2, 7,:= ^i^^ ^ , (42) 



where ujq = v^det K and f^o = V— trK + 2\/det K in the assumption that detK > and trK < 0. Due to 
the singularity the gyroscopic stabilization in the presence of dissipative and non-conservative positional forces 
depends on the ratio and is thus very sensitive to non-Hamiltonian perturbations. We will discuss gyroscopic 
stabilization in more detail in section FZ. II 

We note that the sensitivity of simple eigenvalues of Hamiltonian and gyroscopic systems to dissipative 
perturbations was a subject of intensive investigations, see, e.g., MacKay [BO , Haller [53], and Bloch et al. 
[9]. MacKay pointed out the necessity to extend such a perturbation analysis to multiple eigenvalues [60] . 
Maddocks and Overton [61 initiated the study of multiple eigenvalues and showed that for an appropriate class 
of dissipatively perturbed Hamiltonian systems, the number of unstable modes of the dynamics linearized at a 
nondegenerate equilibrium is determined solely by the index of the equilibrium regarded as a critical point of the 
Hamiltonian. They analyzed the movement of the eigenvalues in the limit of vanishing dissipation without direct 
application, however, to the destabilization paradox and approximation of the singular stability boundary. Our 
calculations performed in this section use the ideas developed in [44l |45l [46l [5^ that, however, can be traced 
back to the works of Andreichikov and Yudovich [IJ and Crandall [19] . 

We see that in Hamiltonian mechanics, the Hamiltonian- Hopf bifurcation in which two pairs of complex 
conjugate eigenvalues approach the imaginary axis symmetrically from the left and right, then merge in double 
purely imaginary eigenvalues and separate along the imaginary axis (or the reverse) has codimension one. In 
the general case of non-Hamiltonian vector fields, the occurrence of double imaginary eigenvalues has codimen- 
sion three. The interface between these two cases possesses the Whitney umbrella singularity; the Hamiltonian 
systems lie on its handle. Quoting Langford from his introductory paper [55] linking Hopf bifurcation, Hamil- 
tonian mechanics and Whitneys umbrella: 'Hopf meets Hamilton under Whitney's umbrella', which, we add, 
was opened by Bottema. 



6 Parametric resonance in systems with dissipation. 

Parametric resonance arises usually in applications if we have an independent (periodic) source of energy. The 
classical example is the mathematical pendulum with oscillating support and a typical equation studied in this 
context is the Mathieu equation: 

X + {uj^ + e cos vt)x = 0. 

In the case of this equation, basic questions are: for what values of the parameters oj, e, v is the trivial solution 
X — X — Q stable or unstable? Another basic question is, what happens on adding damping effects? In the 
theory, certain resonance relations between the frequencies uj and v play a crucial part. See for instance [4], 
[in], [73], [SS] or [H3] and Fig. ^a.) for this classical case. 
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In applications with parametric excitation where usually more degrees of freedom play a part, many combi- 
nation resonances are possible. For a number of interesting cases, analysis and more references see |10l 173] . In 
what follows, the so-called sum resonance will be important. 

First we will consider the general procedure for systems with this combination resonance, after which we 
will discuss an application. 

6.1 Normalization of oscillators in sum resonance 

In |30| a geometrical explanation is presented for damping induced instability in parametric systems using 'all' 
the parameters of the system as unfolding parameters. It will turn out that, using symmetry and normalization, 
four parameters are needed to give a complete description in a two degrees of freedom system, or more generally 
systems where three frequencies are in resonance, but three parameters suffice to visualize the situation. 
Consider the following type of nonlinear differential equation with three frequencies 



which describes for instance a system of two parametrically forced coupled oscillators. A is a 4 x 4 matrix, 
containing a number of parameters, with purely imaginary eigenvalues ±iuji and ±iuj2. Assume that A is 
semi-simple, so, if necessary, we can put A into diagonal form. The vector valued function f contains both 
linear and nonlinear terms and is 27r-periodic in luqI, f(0,a;o^) — for all t. Eq. (|43p can be resonant in many 
different ways, but as announced, we consider here the sum resonance 



where the system may exhibit instability. The parameter S is used to control the detuning d = (61,62) of the 
frequencies {u}i,lu2) near resonance and the parameter fi — (/ii,/i2) derives from the damping coefficients. So 
we may put A = A{6,fi). We summarize the analysis from |30| . 

The basic approach will be to put eq. (I43p into normal form by normalization or averaging whereas the 
theory from [3] will play a part. In the normalized equation the time-dependence is removed from lower order 
and appears only in the higher order terms. It turns out that the autonomous, linear part of this equation 
contains already enough information to determine the stability regions of small amplitude oscillations near the 
origin. The linear part of the normal form can be written as 



X = Ax + ef (x, ujot), X G M^, 



(43) 



Ll)i+UJ2 = Wo, 



z 



A{6, n)z 



with 4-dimensional 




(44) 



and 




(45) 
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Figure 5: The critical surface in (5+) space for eq. (f43|) . /i+ = /ii + /i2, = Mi ~ A*2, (5+ = (^i + (52- 

Only the part fi+ > and (5+ > is shown. The parameters Si, 62 control the detuning of the frequencies, the 
parameters /ii,/i2 the damping of the oscillators (vertical direction). The base of the umbrella lies along the 
(5_i_-axis. 

Since A((5, /i) is the complexification of a real matrix, it commutes with complex conjugation. Furthermore, 
according to the normalization described in [4], [32] and [691 and if uji and uj2 are independent over the integers, 
the normal form of eq. ()43|) has a continuous symmetry group. The second step is then to test the linear part 
A((5, /i) of the normalized equation for structural stability i.e. to answer the question whether there exist open 
sets in parameter space where the dynamics is qualitatively the same. The analysis follows |3] and [4]. The 
family of matrices A((5, fi) is parameterized by the detuning d and the damping fi. The procedure is to identify 
the most degenerate member N of this family, which turns out to be A{d, 0) and then show that A{5, ^) is its 
versal unfolding in the sense of [4j. The family A{S,fi) is equivalent to a versal unfolding of the degenerate 
member N. For details we refer again to [30l[83], an explicit example is discussed in the next subsection. 

We can put the conclusions in a different way: the family A((5, fi) is structurally stable for S, fi > 0, whereas 
A(5, 0) is not. This has interesting consequences in applications as small damping and zero damping may 
exhibit very different behavior. In parameter space, the stability regions of the trivial solution are separated 
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by a critical surface which is the hypersurface where A{S,fi) has at least one pair of purely imaginary complex 
conjugate eigenvalues. As before, this critical surface is diffeomorphic to the Whitney umbrella, see Fig. [5] 
It is the singularity of the Whitney umbrella that causes the discontinuous behavior displayed in the stability 
diagram in the subsection 16.31 The structural stability argument guarantees that the results are 'universally 
valid', i.e. they qualitatively hold for generic systems in sum resonance. 

Above we have described the basic normalization approach, but if we are interested only in the shape of the 
resonance (instability) tongues, there are faster methods. For instance using the Poincare-Linstedt method, see 

m- 

6.2 Rotor dynamics without damping 

The effects of adding linear damping to a parametrically excited system have already been observed and de- 
scribed in for instance [TO], [HS]) [ZH]) or [73]. The following example is based on [55] , 




Figure 6: Rotor with diskmass M, elastically mounted with axial (u) and lateral directions. 

Consider a rigid rotor consisting of a heavy disk of mass M which is rotating with constant rotation speed 
f2 around an axis. The axis of rotation is elastically mounted on a foundation; the connections which are 
holding the rotor in an upright position are also elastic. To describe the position of the rotor we have the axial 
displacement u in the vertical direction (positive upwards), the angle of the axis of rotation with respect to the 
z-axis and around the z-axis. Instead of these two angles we will use the projection of the center of gravity 
motion on the horizontal {x, ?/)-plane, see Fig. [51 Assuming small oscillations in the upright (u) position. 
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frequency 2?], the equations of motion without damping become after rescahng: 

X + 2ay + (1 + Aerj^ cos 2rit)x — 0, 

+ {1 + 467]"^ cos2r]t)y ^ 0. (46) 

The parameter a is proportional to the rotation speed fi. System constitutes a conservative system of 
coupled Mathieu-like equations. Abbreviating P{t) = Arf cos 2rit^ the corresponding Hamiltonian is: 

H = i(l + + eP{t))x' + ^pI + i(l + a^+ eP{t))y^ + + axpy ~ ayp,, 

where Px,Py are the momenta. The natural frequencies of the unperturbed system (j46p . e = 0, are lui = 
\J cP- + \ + a and a;2 = Va^ + 1 — a. By putting z = a; + iy, system (j46p can be written as: 

z - 2aiz + (1 + 4£ry^ cos 2?7i)2 = 0. (47) 

Introducing the new variable: 



e-*"*z, (48) 



and rescaling time rft — r, we obtain: 



4£cos2r U = 0, (49) 



where the prime denotes differentiation with respect to t. By writing down the real and imaginary parts of this 
equation, we have actually got two identical Mathieu equations. 

Using the classical and well-known results on the Mathieu equation, we conclude that the trivial solution is 



stable for e small enough, provided that vT+c? is not close to nrj , for n = 1, 2, 3, .... The first-order and most 
prominent interval of instability, n — 1, arises if: 



\/l-t-a2 « ry. (50) 

If condition (l50|) is satisfied, the trivial solution of equation ((49|) is unstable. Therefore, the trivial solution of 
system (|46p is also unstable. Note that this instability arises when: 

uji+u)2 = 2r], 

i.e. when the sum of the eigenfrequencies of the unperturbed system equals the excitation frequency 2rj which 
is the sum resonance of first order. The domain of instability is bounded by: 



VT+^ {I ± e) + 0{e^) . (51) 



See Fig. [7]Jb) where the V-shaped instability domain is presented in the case of rotor rotation {a ^ 0) without 
damping. 

Higher order combination resonances can be studied in the same way; the domains of instability in parameter 
space continue to narrow as n increases. As noted, the parameter a is proportional to the rotation speed f2 of 
the disk and also to the ratio of the moments of inertia. 
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6.3 Rotor dynamics with damping 

We add small linear damping to system (j46p . with positive damping parameter fj, = 2eK. This leads to the 
equations: 

X + 2ay + (1 + 4e77^ cos 2r]t)x + 2eKx = 0, 

y - 2ax + (1 + Aeri'^ cos 2r]t)y + 2eKy = 0. (52) 
and using the complex variable z: 

z — 2aiz + {l + 4:eT]'^cos2r]t) z + 2eKz^0. (53) 

Because of the damping term, we can no longer reduce the complex eq. (j53p to two identical second order 
real equations, as we did previously. 

In the sum resonance of the first order, we have uji + lu2 ~ 2r] and the solution of the unperturbed (e = 0) 
equation can be written as: 

zit) = zie"^^' + Z2e-'"^*, zi, Z2 e C, (54) 



with uji = \/(?+T + a, L02 — y oP' + 1 — a. 

Applying variation of constants leads to equations for zi and Z2'. 

ie 



zi = 



Z2 = 



-{2K{iuJiZi - iuj2Z2e-'^'^^+'^^^^) + 



t^l + UJ2 

Aif cos2ryt(zi + Z2e"'('^i+""^*)), 



-{2K{iuJlZie'^'^'+'^^'>' - iuj2Z2) + 



UJi + U}2 

4772 cos 2r/i(zie*("'^+'^^)* +22))- (55) 

To calculate the instability interval around the value 770 = hi^i + ^2) = vc? + l, we apply normal form or 
(periodic solution) perturbation theory, see [SS] for details, to find for the stability boundary: 



rib = y/TTa^ il±sJl + a^ - ^ 



= Vl + a^ yi±^{l + a2)£2 _ ^ j (56) 

It follows that, as in other examples we have seen, the domain of instability actually becomes larger when 
damping is introduced. See Fig. [TId. 

The instability interval, shows a discontinuity at k = 0. 

If K — > 0, then the boundaries of the instability domain tend to the limits t]i, y/1 + a'^{l ± evT+7?) 



which differs from the result we found when k — : 7/h = vT+c?(l ± e). For reasons of comparison, we display 
the instability tongues in Fig. [7] in the four cases with and without rotation, with and without damping. 

Mathematically, the bifurcational behavior is again described by the Whitney umbrella as indicated in 
subsection 16. II In mechanical terms, the broadening of the instability-domain is caused by the coupling between 
the two degrees of freedom of the rotor in lateral directions which arises in the presence of damping. 
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Figure 7: (a) The classical case as we find for instance for the Mathieu equation with and without damping; in 
the case of damping the instability tongue is lifted off from the ?7-axis and the instability domain is reduced, 
(b) The instability tongues for the rotor system. Again, because of damping the instability tongue is lifted off 
from the 77-axis, but the tongue broadens. The boundaries of the ^-shaped tongue without damping are to first 
approximation described by the expression 77 = VT + c?(l ± e), 770 — vl + Q?. 

7 Manifestation of the destabilization paradox in other apphcations 

In this section we discuss additional applications from physics and engineering, both finite- and infinite- 
dimensional. 

7.1 Gyroscopic systems of rotor dynamics 

Investigation of the stability of equilibria of the Hanger's [21] and Crandall's [HI [67] gyropendulums as well as 
of the Tippe Top [TH [S3] and the Rising Egg [T5] leads to the system of linear equations known as the modified 
Maxwell-Bloch equations [9]. 

The modified Maxwell-Bloch equations are the normal form for rotationally symmetric, planar dynamical 
systems [^[[TS]. They follow from the equation ([25)1 for m = 2, D = I, and K = kI, where k, corresponds to 
potential forces, and thus can be written as a single differential equation with complex coefficients 

X -\- ifii + (Ji + ivx + Ka; = 0, x = x\ — 1x2- (57) 

According to ([TTj) the solution .t = of equation (|57|) is asymptotically stable if and only if 

(5 > 0, ^>---K. (58) 

V 

For K > the domain of asymptotic stability is a dihedral angle with the fi-axis serving as its edge, Fig.[8l^b). 
Its sections by the planes il = const are contained in the angle-shaped regions with the boundaries 

V ^ ^ ^- (59) 
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At Q = the angle is bounded by the hnes v — ±(5y^ and thus it is less than tt. The domain of asymptotic 
stability is twisting around the fJ-axis in such a manner that it always remains in the half-space (5 > 0, Fig.[5I^b). 
Consequently, the system that is statically stable at = and (5 > can become unstable at greater f2 in the 
presence of non-conservative positional forces, as shown in Fig. [U^b) by the dashed line. The larger magnitudes 
of circulatory forces, the lower |f2| at the onset of instability. This is a typical example of dissipation-induced 
instability in the sense of [HI 1521 1531 154] when only non-Hamiltonian perturbations can cause the destabilizing 
movements of eigenvalues with definite Krein signature [48] . 

T\ an 



Figure 8: (a) Hanger's gyropendulum; (b) Dissipation-induced destabilization of its statically stable equilibrium 
(k > 0) in the presence of circulatory forces; (c) Singular domain of gyroscopic stabilization of its statically un- 
stable equilibrium (k < 0) in the presence of non-Hamiltonian perturbations yields the destabilization paradox. 

As K > decreases, the hypersurfaces forming the dihedral angle approach each other so that, at k = 0, they 
temporarily merge along the line — and a new configuration originates for k < 0, Fig.[5]Jc). The new domain 
of asymptotic stability consists of two disjoint parts that are pockets of two Whitney's umbrellas singled out by 
inequality S > 0. The absolute values of the gyroscopic parameter Q in the stability domain are always not less 
than f^o = 2\/—K. As a consequence, the system that is statically unstable at f2 = can become asymptotically 
stable at greater il in the presence of circulatory forces, as shown in Fig. [8t^c) by the dashed line. 

As a mechanical example we consider Hanger's gyropendulum [241 , which is an axisymmetric rigid body 
of mass m hinged at the point O on the axis of symmetry as shown in Fig. ©(a). The body's moment of 
inertia with respect to the axis through the point O perpendicular to the axis of symmetry is denoted by /, 
the body's moment of inertia with respect to the axis of symmetry is denoted by Iq, and the distance between 
the fastening point and the center of mass is s. The orientation of the pendulum, which is associated with 
the trihedron OxfUfZf, with respect to the fixed trihedron OxiUiZi is specified by the angles ?/;, 0, and (j). The 
pendulum experiences the force of gravity G — mg and a follower torque T that lies in the plane of the Zi and 
Zf coordinate axes. The moment vector makes an angle of 770; with the axis Zi, where 77 is a parameter (?7 7^ 1) 
and a is the angle between the Zi and z/ axes. Additionally, the pendulum experiences the restoring elastic 
moment R — ^ra in the hinge and the dissipative moments B — —bus and K — ~k(j), where ujs is the angular 
velocity of an auxiliary coordinate system OxsVsZg with respect to the inertial system and r, 6, and k are the 
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corresponding coefficients. 

Linearization of the nonlinear equations of motion derived in [24j with the new variables xi — ip and X2 = d 
and the subsequent nondimensionalization yield the Maxwell-Bloch equations (j57p where the dimensionless 
parameters are given by 



The domain of asymptotic stability of the Hauger gyropendulum, given by (j58p . is shown in Fig. [5fb,c). 

For the statically unstable gyropendulum (k < 0) the singular points on the fi-axis correspond to the critical 
values ±r2o = ±2^/— k and the critical frequency o^o = V"— k. We find approximations of the stability boundary 
near the Whitney umbrella singularity as derived in [46|, I50j : 



n,r{i^,s) = ±2V^± (61) 



Thus, Hanger's gyropendulum, which is statically unstable at il = 0, can become asymptotically stable for 
sufficiently large > Qq under a suitable distribution of dissipative and nonconservative positional forces. 
For almost all combinations of S and v the onset of gyroscopic stabilization of the non-conservative system is 
greater than that of a pure gyroscopic one (destabilization paradox: Qcrij^iS) >^o)- The obtained results are 
vahd also for the equihbria of Tippe Top, Rising Egg, and Crandall's gyropendulum [44),l45]. 



7.2 Circulatory systems of rotor dynamics 

In some rotor dynamics applications gyroscopic effects are neglected |34 [ I53 j . For example, in the modeling 
of friction-induced oscillations in disc- and drum brakes, clutches and other machinery, the speed of rotation 
is assumed to be small. This frequently yields the linearized equations of motion in the form of a circulatory 
system with or without damping. In recent models the damping is included because it is believed that high 
sensitivity of the squeal onset to the damping distribution might be responsible for the poor reproducibility of 
the laboratory experiments with the squealing machinery. 

Hoffmann and Gaul 28J studied a model of a mass sliding over a conveyor belt with friction and detected 
that small damping in this circulatory system destroys the reversible Hopf bifurcation and makes the collision 
of eigenvalues imperfect, exactly as it happens with the eigenvalues of Ziegler's pendulum |401 147j . 

In order to study squeal vibration in drum brakes Hulten [31[ I75j introduced a model shown in Fig. [9ja). 
This model is composed of a mass m held against a moving band; the contact between the mass and the band 
is modeled by two plates supported by two different springs. It is assumed that the mass and band surfaces are 
always in contact and that the contact can be expressed by two cubic stiffnesses. Damping is included as shown 
in Fig. [IJa). The friction coefficient at contact is assumed to be constant and the band moves at a constant 
velocity. Then it is assumed that the direction of friction force does not change because the relative velocity 
between the band speed and ii or ±2 is assumed to be positive. The tangential force Ft due to friction contact 
is assumed to be proportional to the normal force as given by Coulomb's law: Ft = ^J.Fn- Assuming the 
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Figure 9: (a) A model of a drum brake [3T]; (b) Its numerically calculated critical friction coefficient at the 
onset of flutter instability as a function of damping parameters [75] ; (c) A model of a disc brake (QA^ ; (d) Its 
critical friction coefficient at the onset of flutter instability as a function of damping parameters calculated both 
numerically and perturbatively [47j . 



normal force Fpf is linearly related to the displacement of the mass normal to the contact surface, the resulting 
equations of motion can be expressed as 










: 













X = 0, 



?72W0,2 / \ ^J■U!o^l 



■^0,2 



being exactly of the form considered by Bottema. Here the relative damping coefficients are denoted by rji = 
Ci/y/rriiki {i — 1,2) and natural pulsations are wo.i = \/ /™i (* — lj2). Fig. [Sl^b) shows the numerically 
calculated domain of asymptotic stability of the drum brake in the space of the friction coefficient /x and two 
damping coefficients ryi and r]i with the Whitney umbrella singularity [75j . 

In Fig. in^c) a model of a disc brake proposed in [64] is demonstrated. Its linearized equations of motion are 
again that of a circulatory system with small damping. It is not surprising that the critical friction coefficient 
at the onset of friction-induced vibrations as a function of two damping coefficients is represented in Fig. [9l[d) 



24 



by a surface with the Whitney umbreUa singularity [47] . 

In both examples a selected distribution of damping exists that yields an increase in the critical load rather 
than decrease that happens for all other distributions. This possibility for stabilization was first pointed out in 
[71] for the Ziegler's pendulum. We will discuss this effect below in more detail. 

7.3 Infinite-dimensional near-reversible and near-Hamiltonian systems 

Dynamic instability, or flutter, is a general phenomenon which commonly occurs in coupled fluid-structure 
systems including pipes conveying fluids and airfoils [TOl [26l [27] . Typically, the models are finite dimensional 
or continuous reversible systems that demonstrate the destabilization paradox in the presence of damping. In 
a recent study ^87j Ziegler's paradox was observed in a problem of a vocal fold vibration (phonation) onset. 

7.3.1 Near-reversible case: Beck's column with external and internal damping 

Beck's column loaded by a follower force is a paradigmatic model for studying dynamical instability of struc- 
tures. In 1969 Bolotin and Zhinzher llj investigated the effects of damping distribution on its stability. They 
considered on the interval x G [0, 1] the non-selfadjoint boundary eigenvalue problem of the form 

Lu := N{q)u + AD(rfi, da)" + A^Mu = 0, Uu := UAr(9)u + AU£,(di, d2)u + A^U^u = 0, (62) 

where A is an eigenvalue with the eigenfunction u{x). The operators in the differential expression 

N = dt + qdl, D = did^ + d2l, M = I (63) 

depend on the magnitude of the follower load q and the parameters of external, rfa, and internal (Kelvin- Voight), 
di, damping. The matrices of boundary conditions in |11J are \] d = 0, Um = 0, and 

/lOOOOOOO^ 
1 
1 

yoooooooiy 

and the vector u = (m(0), 9a;u(0), 9^u(0), 9^u(0), m(1), 9a;M(l), 9^u(l), 9^u(l))-^. Some authors considered differ- 
ent boundary conditions that depend both on the physical parameters and on the spectral parameter [63[ 188] 

The undamped Beck's column is stable for q < qo ~ 20.05 [T^. Stability is lost at g > go when after 
the reversible Hopf bifurcation the double pure imaginary eigenvalue iwo — 11.02 splits into a pair of complex 
eigenvalues. In jll, it was found that in the presence of infinitesimally small Kelvin- Voight damping the critical 
load is reduced to q = q^r — 10.94 and the critical frequency drops to uj = ojcr — 5.40. 

There were numerous attempts to find an approximation of the new critical load by studying the splitting of 
the double eigenvalue zwq of the unperturbed reversible system due to dissipative perturbations [70] . Banichuk 
et al. [H] [7j have emphasized the importance of degenerate perturbations, the linear part of which is in the 
tangent plane to the Whitney umbrella singularity. Nevertheless, their analysis is not complete. 



(64) 
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Further development of the approach of [51 [7] in [371 HOI HIl HH 1131 IZl] resulted in the approximation to the 
critical load in the form 

/(h,d)2 / 

where the vector of the damping parameters d=((ii, 1^2) and angular brackets denote the scalar product in K.^. 
The components of the vector f and the real scalar / are 

and the components of the vector h are defined as 

with the asterisk denoting complex conjugate transposition and (w, w) = u{x)v{x)dx. The derivatives are 
taken at d = and q — qq corresponding to the eigenvalue A — iluq with the eigen- and associated functions uq 
and ui. The real matrix H has the components 

and the real matrix G is defined by the expression 

(«d^d)=t'^^(|£-2'-o+vSVS^^.), (69) 

r— 1 ^ ^ 

where W2 is the solution of the boundary value problem 

N{qo)w2 - UJ0MW2 = ^ ^''a^^"' U^(^o)^2 - ^oUAfW2 = ^ 

r— 1 ^ r— 1 ^ 

The eigenfunctions uq and wq and the associated functions ui and wi of the original and adjoint eigenvalue 
problems are chosen to satisfy the bi-orthogonality and normalization conditions 

2iujo{Mui,vi) + {MuQ,vi) + {Mui,vq) + (VqVi + V^vo)*(2itJoUAfUi + UmUq) + v*V*Ua/Ui = 0, 
2iLUo{Mui,vo) + 2it^oVoV*UMUi = 1, (71) 

where the adjoint boundary value problems are connected by the Lagrange formula 

(Lu, v) - {u, L*v) = (Vv)*Uu - (Vv)*Uu. (72) 

Formula (|65p can serve for the approximation of the jump in the critical load. In the finite dimensional case 
with two degrees of freedom the expression for the limit of the critical load reduces to p4|) [40] , For the Beck 
column described by the equations (|62|) we calculate the critical load as [HI |42] 



1 009(72 

'^"^'^^'^^ - '^^ - (14.34d, + 0.091ci.)2 + '^-'^"^'^ + (^3) 
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Additionally, (j65p captures more information on the geometry of the stability domain. For example one can plot 
the cross sections of the stability domain ([75]) for the different levels of q and find that some combinations of 
internal and external damping increase the critical load for the Beck's column, as shown in the central and right 
pictures of Fig. 1101 The form of the stability boundary with the Whitney umbrella singularity approximated 
by equation (|73p was confirmed later by numerical computations in [33 . The limit in the critical load following 
from ([75)1 agrees well with the numerical data of [T] . 



0.2 

qcr= qo - 5 



0.4 





Figure 10: (Hatched) Cross sections of the approximation ([75| to the stability domain of the damped Beck's 
column showing the possibility of the increase of the critical load by small damping. 



Structural mechanics also has examples of near-Hamiltonian continuous systems showing discontinuous 
changes in the stability domain. As a modern application we mention a moving beam with frictional con- 
tact investigated in [77| . Below we will consider an interesting example of the occurrence of the destabilization 
paradox in fluid dynamics. 



7.3.2 Near-Hamiltonian case: The instability of baroclinic zonal currents 

In the 1940s the first studies appeared of instability of baroclinic zonal (west-east) currents in the Earth's 
atmosphere [HI [20] . It is remarkable that the unexpected destabilizing effect due to the introduction of dissi- 
pation was discovered in the linear stability analyses of this hydrodynamical problem by Holopainen (1961) |29j 
and Romea (1977) [65, at the very same period of active research on the destabilization paradox in structural 
mechanics. Recently these studies were revisited by Krechetnikov and Marsden [54j with the aim to handle 
rigorously the treatment of dissipation-induced instability. 

Romea considered an infinite channel in the periodic zonal direction x of width L in the meridional direction 
y that is rotating with an angular velocity Q. Two layers of incompressible, homogeneous fluids of slightly 
different densities (the lighter fluid on top) are confined by the side walls and by horizontal planes, a distance D 
apart. For simplicity, it is assumed that, in the absence of motion, the interface is located halfway between the 
horizontal planes, and is fiat so that centrifugal effects may be ignored. Each layer moves downstream with a 
constant velocity and the slope of the interface is related to these velocities through the thermal wind relation. 
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It is implicitly assumed that this basic state is maintained against dissipation by an external energy source 
which is unimportant with respect to the rest of the problem [65) . 

The linearized equations for each layer near the basic state, characterized by the geostrophic streamfunctions 
~Uiy and ~U2y, are according to [65ll54] : 



{dt + UA)[W^^i+F{ip2-^i)] + [fi + F{Ui-U2)]d^V^i = -rVVi, 

{dt + U2d^)[V^ip2 + F{ipi-ip2)] + [P-FiUi~U2)]d,iP2 = -rV^iP2- (74) 

where F is the internal rotational Froude number, r > is the measure of the effect of Ekman suction (Ekman 
layer dissipation), and /? is the planetary vorticity factor introduced to take into account the variation of the 
Coriolis parameter with latitude (/3-effect). 

Assuming the wave solutions ipi^2 ~ e*"^^"*^'-* sin(m7r?;), where real a > is the x wavenumber, Romea 
obtained a dispersion relation for the complex phase speed c = + ict in the form of the second-order complex 
polynomial. The real part of c is the speed of propagation of the perturbation, while aci is the growth rate of 
the wave. If Ci > 0, the wave grows, and the system is unstable. 

In the inviscid case when the Ekman layer dissipation is set to zero, the transition to instability occurs 
through the Krein collision that occurs at Uc ■— Ui — U2 = Ud, where [551 [Mj 

t/e. = (75) 

with = + to^tt^. The critical shear Ud as a function of the wavenumber is plotted in Fig. [TlTleft). This 
curve bounds the region of marginal stability of the system without dissipation. 

In the limit of vanishing viscosity (r ^ 0), the stability boundary differs from (jTS]) 

U^^ = (76) 

a(a2 + F)V2F - ' 

The discrepancy between the stability domains of viscous and inviscid systems is clearly seen in Fig. [TTf left). 

Therefore, Romea demonstrated that an introduction of infinitesimally small dissipation destabilizes the system, 

lowering the curve of marginal stability by an 0(1) amount. This is the appearance of the destabilization paradox 

in a continuous near-Hamiltonian system, which is similar to that found in near-reversible systems like Ziegler's 

pendulum (cf. Fig. [3]) and Beck's column with dissipation [TTl[80]. Fig. fTlT right) reproduces the original drawing 

from [65 showing the typical imperfect merging of modes that substitutes the 'perfect' Krein collision in 

near-Hamiltonian and near-reversible systems. Approximation to the eigenvalue branches in imperfect merging 

can be efhciently calculated by means of the perturbation theory of multiple eigenvalues for a wide class of 

non-conservative systems [37 l [39 l l40 l l42 l [43l [49] . 

8 Conclusion 

We have revisited the pioneering result of Oene Bottema who in 1956 resolved the paradox of destabilization by 
small damping and interpreted it by means of what is now called the Whitney umbrella singularity. We have 
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0.5 1.0 1.5 2.0 

Figure 11: Original drawings from the 1977 work of Romea [65]: (left) Critical shear as a function of wavenum- 
ber demonstrates a discontinuous transition from the case when the Ekman layer dissipation r — initially 
(Uci) to the case when r ^ (Ucr)', (right) a typical imperfect merging of modes (growth rates) that substi- 
tutes the 'perfect' Krein collision in near-Hamiltonian and near-reversible systems and is characteristic for the 
destabilization paradox [TT], [T^ . 

shown that this phenomenon frequently occurs in near-Hamiltonian and near-reversible systems originating in 
very different areas of mechanics and physics ranging from hydrodynamics to contact mechanics and we have 
presented a unified treatment of it. There are a few related topics upon we did not touch. We mention interesting 
connections of this effect to structured pseudospectra j35] and to eigenvalue optimization problems [1^. We 
did not even consider the effect of nonlinearites. We mention the closely related effect of discontinuous change 
of the critical flutter frequency due to small dissipation [11] |42] and its connection to the Whitney umbrella 
singularity at the exceptional points on the eigenvalue surfaces [75] • Another related topic is the role of the 
spectral exceptional points in modern non-Hermitian physics including crystal optics, open quantum systems, 
and T'T-symmetric quantum mechanics [8 . All this shows that modern non-conservative and non-Hermitian 
problems are a perfect field of applied mathematics with a big potential for new discoveries. 
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